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ABSTRACT 

We explore the properties of high-redshift Lyman-alpha emitters (LAE), and 
their link with the Lyman-Break galaxy population (LBG), using a semi-analytic 
model of galaxy formation that takes into account resonant scattering of 
Lya photons in gas outflows. We can reasonably reproduce the abundances of 
LAEs and LBGs from z « 3 to 7, as well as most UV LFs of LAEs. The stronger 
dust attenuation for (resonant) Lya photons compared to UV continuum photons 
in bright LBGs provides a natural interpretation to the increase of the LAE frac¬ 
tion in LBG samples, Alae, towards fainter magnitudes. The redshift evolution 
of Alae seems however very sensitive to UV magnitudes limits and EW cuts. In 
spite of the apparent good match between the statistical properties predicted by 
the model and the observations, we find that the tail of the Lya equivalent width 
distribution (EW > 100 A) cannot be explained by our model, and we need to 
invoke additional mechanisms. We find that LAEs and LBGs span a very similar 
dynamical range, but bright LAEs are ~ 4 times rarer than LBGs in massive 
halos. Moreover, massive halos mainly contain weak LAEs in our model, which 
might introduce a bias towards low-mass halos in surveys which select sources 
with high EW cuts. Overall, our results are consistent with the idea that LAEs 
and LBGs make a very similar galaxy population. Their apparent differences 
seem mainly due to EW selections, UV detection limits, and a decreasing Lya to 
UV escape fraction ratio in high SFR galaxies. 


Key words: galaxies: formation - galaxies: evolution - galaxies: high-redshift - 
radiative transfer - methods: numerical. 


1 INTRODUCTION 


Star-forming galaxies are commonly detected nowadays 
via two main channels at high redshift (z > 3). On the 
one hand, the Lyman-Break technique is quite efficient at 
selecting objects with strong stellar UV continuum (A « 
1500 A), using a set of broad band colour-colour criteria 
(e.g. Steidel et al.||1996 1999 2003 Gabasch et al.|[2004| 


[Bouwens et al.||2007 20111. On the other hand, narrow- 

band imaging and spectroscopic surveys have also been 
able to probe large samples of high redshift sources via 
their (nebular) Lya emission line (A = 1215.67 A; e.g. 


Hu et al.|1998||Rhoads et al.|2000||Shimasaku et al. 2006| 

Ouchi et al. 

20081 Rauch et al.|20081 Hu et al.|20101|Hayes| 

et al.|2010| 

Gassata et al.|2011| Blanc et al.|2011|). 
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The joint study of these two populations is essen¬ 
tial to improve our understanding of the formation and 
evolution of galaxies in the early Universe, and their im¬ 
pact on the intergalactic gas notably during the epoch of 
reionisation. Large observational datasets have enabled 
to constrain many statistical properties of Lyman-Break 
Galaxies (LBG) and Lyman-alpha Emitters (LAEs). The 
evolution of the UV luminosity function (LF) of dropout 
galaxies appears to be much stronger than the Lya LF of 
LAEs at z ~ 3-7. In terms of physical properties, typ¬ 
ical LBGs correspond to massive objects, with stellar 


mass correlated to UV luminosity (Shapley et al. 2001 


[Bouwens et al.|2010||Gonzalez et al.|2011 1. LAEs are com¬ 
monly thought to be less massive galaxies and to form a 
highly inhomogeneous population in terms of age, mass 


or dust content ( 

Gawiser et al. 

2006 

Finkelstein et al. 

2007 

IPirzkal et al.|2007||Finkelstein et al.|2009t|0no et al. 


20101. Although part of the difference between the LBG 
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and LAE populations is inherent to their respective, both 
biased, methods of selection, their link is not quite well 
established yet, and it is still unclear to what extent they 
are representative of the underlying galaxy population. 

In the modern picture of hierarchical structure forma¬ 
tion, galaxies form through gas accretion within virialised 
halos of dark matter located at the density peaks of the 
background matter field | Birnboim k, Detel]|2003[ |Keres| 


galaxies ar 
halos. The 
to be one 
rently obse 
tion measi 

e tightly linked to the characteristics of the host 
host halo masses of typical LBGs are inferred 
order of magnitude larger than those of cur- 
rved LAEs from angular auto-correlation func- 

irements (Hamana et al. 2004 Gawiser et al. 

2007| [Hildebrandt et al.|2009| [McLure et al.[[2009| [Ouchi 

et al. 

2010 

and halo occupation distribution models (Lee 

et al. 

2006 

[Jose et al.|2013[). 


Accurate estimates on the halo masses of LBGs and 
LAEs at high-redshift are crucial to constrain formation 
history of local galaxies and identify the progenitors of 
the Milky-Way in the hierarchical context. The observed 
LBGs and LAEs luminosities span several orders of mag¬ 
nitude and their host halo properties are thus expected to 
cover a wide range. Although it is essential to investigate 
the connection between UV/Lya-selected galaxies and 
their dark matter halos within cosmological simulations, 
special care has first to be taken to describe the complex 
transfer of Lya photons in the interstellar medium. 

The line emission originates in Hll regions in the in¬ 
terstellar medium as ionising radiation produced by mas¬ 
sive, short-lived, stars is reprocessed into Lya photons. 
The intrinsic Lya emission is then a direct tracer of re¬ 
cent star formation activity in galaxies. However, the ob¬ 
served flux can be strongly reduced by dust extinction in 
the interstellar medium (ISM) and intervening hydrogen 
absorption in the intergalactic medium (IGM). Due to a 
large absorption cross-section combined with substantial 
Hi column densities, the medium becomes optically thick 
to Lya photons which undergo a resonant scattering pro¬ 
cess. Lya photons diffuse both in frequency and physical 
space which can strongly alter the shape of the observed 
line. The travelling path in neutral gas, as well as the 
probability of encounter with dust grains, may be dra¬ 
matically increased, and interactions with atoms in the 
tail of the velocity distribution can scatter photons off to 
the wing. Lya radiation transfer is thus highly sensitive 
to the geometry, the ionisation state, and the kinematics 
of the ISM. 

Only very idealised cases can be investigated ana¬ 


lytically (e.g. Neufeld 1990 Dijkstra et al. 20061 and 
more realistic configurations require numerical schemes. 


mostly based on the Monte-Carlo technique ( 

Ahn et al. 

2001 

[Zheng & Miralda-Escude 2002| [Verhamme et al. 

2006 

[Dijkstra et al. 2006| Hansen & Ohj 2006| Laursen 

& Sommer-Larsen 2007[ Laursen et al. |2013|). Although 


a plethora of codes have been developed, only a few au¬ 
thors have intended to address the Lya RT issue on galac¬ 


tic scales within hydrodynamical simulations (Tasitsiomi 
2006| [Laursen et al.|[^09[ IBar nes et aL]|2011| |Yajima' 
et al.|[^12| [Verhamme et al. [2012 1. These studies are ex¬ 
tremely useful because they can follow the propagation of 
Lya photons through more realistic, non-linear, ISM den¬ 
sity/velocity fields. Even though some trends are identi¬ 
fied (anti-correlation between halo mass and Lya escape 
fraction, orientation effects, etc), the results still appear 
to depend on the resolution and the physics implemented 


in the underlying hydrodynamical simulation (e.g. [VeF^ 


hamme et al. 20121. Moreover, only a handful of high- 


resolution galaxies can be studied at once due to comput¬ 
ing time issues. 

Semi-analytic models (SAM) of galaxy formation 
provide a unique alternative to investigate statistical sam¬ 
ples of mock galaxies and their evolution with redshift. 
Recently, several models accounting for the Lya RT in 
galaxies based on numerical libraries have been presented. 
Rather simplistic physical pictures are usually adopted. 


2012 


such as plane-parallel slab geometries | Forero-Romero 
et al.| 2011 \ or thin expanding shell models (|Qrsi et al. 


Garel et al.||201^. Nonetheless, these approaches 


provide a more accurate treatment than previous mod¬ 
els which used simple parametrisations to describe the 
Lya transfer effects and the observed Lya properties of 
galaxies, assuming a constant Lya escape fraction (e.g. 
|Le Delliou et al.|2005[[Nagamine et al.|2010[ ), simple dust 
models neglecting Lya resonant scattering (|Mao et al.| 


2007 \, or using phenomenological prescriptions to account 
for its effect ||Kobayashi et al.||2007| pSTo) |Dayal et al.| 

[MT| . 

In the present paper, we use the semi-analytic model 
presented in [Garel et al.| ( |2012| ) to investigate the UV, 
Lya, and dark matter halos properties of high redshift 


galaxies. In Garel et al. (2012 hereafter Paper I), we de 


scribed the coupling of GALICS {GALaxies In Cosmo¬ 
logical Simulations; IHsitton et S)]|2003| , a hybrid model 
of galaxy formation based on a N-body cosmological sim¬ 
ulation, with the library of [Schaerer et ah] ( [2011[ ) which 
computes numerically the Lya RT through expanding gas 
shells. We present an overview of our model in Section 
In Section we show that this model is able to repro¬ 
duce the observed UV and Lya luminosity functions from 
z ~ 3 to z ~ 7. In Section and we investigate the 
observed cross properties of LAEs and LBGs to study the 
connection between these two populations, focussing on 
the Lya equivalent width distributions. In Section we 
present our model predictions for the host halo properties 
of LAEs and LBGs in terms of halo mass and occupa¬ 
tion number. Finally, Section [T] gives a discussion and a 
summary of our results. 

All quantities used throughout this paper assume 
the following cosmological parameter values: h = 
i7o/(100kms“^ Mpc"^) = 0.70, Qa = 0 .72, = 0.28, 

fib = 0.046, and erg = 0.82 (WMAP-5; Komatsu et al. 


20091. All magnitudes are expressed in the AB system. 


2 MODEL OVERVIEW 

Our modelling of the high-redshift galaxies and their 
Lya emission properties has been fully described in [Pa-[ 
[per I[ Here, we briefly summarise the main aspects of our 
approach. 

We use an updated version of the GALICS semi- 
analytic galaxy formation model to predict the abun¬ 
dances and physical properties of galaxies, based on the 
hierarchical evolution of dark matter halos. The forma¬ 
tion and evolution of dark matter structures is described 
by a cosmological N-body simulation post-processed with 
a Friends-of-Friends algorithm ( Davis et al.|1985 1 to iden¬ 
tify bound regions (i.e. dark matter halos). The merging 
histories of dark matter halos are computed according to 


the procedure of Tweed et al. (20091. Our N-body simu¬ 
lation has been run with GADGET ( Springel|[2005 I using 
1024® particles in a periodic, comoving box of 100/i“® Mpe 
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on a side. The mass resolution (Mhaio ^ 2 x 10® Mq , cor¬ 
responding to bound groups of at least 20 particles) and 
the box size have been chosen to allow us to investigate 
the statistical properties of LAEs and LBGs currently de¬ 
tectable at high-redshift. While in |Paper~I| we assumed pa¬ 
rameter values from the WMAP-3 release ( |Spergel et al.| 
|2007[ ), here we adopt a more recent set of cosmological 
parameters consistent with the WMAP-5 data, which are 
described above. 

Our modelling of galaxies and Lya emission is very 
similar to Paper I, except for a few changes that we detail 
below. In [Paper l[ we presented the modifications made 
to the original GALICS model described in [Hatton et al.j 
(20031 (see also Cattaneo et al.|2006 l. The main updates 
consisted of a new implementation of gas accretion onto 
galaxies (cold flows), dust extinction of the UV contin¬ 
uum (spherical geometry so as to be consistent with the 
Lya radiative transfer within outflowing shells of gas - see 
below), and star formation. We compute star formation 
rates directly from the cold gas surface density using the 
Kennicutt-Schmidt law (i.e. Esfr oc Eg^^), and we had to 
increase the star formation efficiency by a factor of 25 in 
order to match the observed high redshift UV luminos¬ 
ity functions. This departure from the z=0 normalization 
of |Kennicutt| ( |1998P was certainly due to the underlying 
(WMAP-3) cosmology that we adopted for the dark mat¬ 
ter simulation, and especially the low value of erg (0.76; 
[Spergel et al.|[^07| ) which delays the growth of structure 
at early times. To this extent, boosting star formation 
was necessary to compensate the low abundance of star 
forming galaxies at high-redshift. With our new N-body 
simulation, ng is slightly higher which eases the formation 
of larger structures at earlier epochs. We find that an in¬ 
creasing factor for a star formation efhciency equal to 5 
in our model is sufficient to obtain a reasonable match to 
the high redshift data (see section]^. 

We use the stellar libraries of IDevriendt et ah] 
( 1999[ STARDUST) to compute the intrinsic UV and 
Lya emission of galaxies for a Kennicutt IMF (0.1- 
120 Mq ; |Kennicutt||1983|). Und er the case B approxi¬ 
mation (Osterbrock & Ferland 20061, two third of the 


ionising photons emitted by stars are reprocessed into 
Lya photons when HIT regions recombine, and we assume 
that no ionising photon can escape the medium (optically 
thick limit).The intrinsic Lya luminosity of galaxies is 
therefore given by = 0.67Q{H)^^, where Q{H) is 
the production rate of hydrogen-ionising photons, c the 
speed of light, hp the Planck constant, and Aq, is the Lya 
wavelength at line center. For the shape of the intrinsic 
line, we assume a Gaussian profile centered on Aa with a 
characteristic width given by the circular velocity of the 
galaxy disc ( Santos|2004 1. 


We take into account the transfer of UV continuum 
and Lya radiation through outflowing gas and dust. This 
picture is motivated by the apparent ubiquity of galac¬ 
tic winds seen at high redshift (powered by supernovae 
and/or stellar winds), as well as their impact on the ob¬ 
served Lya spectral signature in Lyman Break galaxies 


and Lyman Alpha Emitters ([Shapley et al. 2003| 

Steidel 

et al.|2010| 

McLinden et al.|2011| 

Finkelstein et a 

l.|2011[ 

Kulas et al. 

2012 

Berry et al.|2012 

1. Following the work of 

Verhamme et al. 

(20061, we describe galactic outflows as 


homogenous, spherical, expanding shells of cold gas mixed 


with dust. Verhamme et al. (20081 have shown that high 


redshift Lya line profiles can be well reproduced by 3D 
Monte-Carlo radiation transfer (RT) models in expand¬ 


ing shells, when adjusting the values of various relevant 
shell parameters: the bulk motion expansion velocity 14xp, 
the gas velocity dispersion h, the dust opacity Tdust, and 
the neutral hydrogen column density of the shell Nm. We 
compute shell parameters for GALICS galaxies in a simi¬ 
lar fashion to [PapefT] In our model, the shell speed V/xp 
has a (weak) dependence on the star formation rate, the 
value of b is fixed to 20 km s”'^, and the dust opacity is 
estimated from the [Mathis et al.| ( |1983[ ) extinction curve 
at a given wavelength, the gas metallicity and column 
density, and a redshift-dependent dust-to-gas ratio (see 
[Garel et al.|201^ for a detailed description of our model). 
The column density is proportional to Mgj^°jj/(47rRsheii®), 
where M®j(gjj and Rsheii are the shell mass and size re¬ 
spectively. We use the cold gas in galaxies as a proxy for 
the shell mass, and the disk scalelength Rdisk to compute 
Rsheii such that Rsheii = CRdisk where ^ is a free param¬ 


eter of the order of unity. While ^ was 1.0 in [Paper l| 
we now set the value to 1.2 to improve the agreement of 
our model with the luminosity functions data over our 
redshift range of interest. 

The escape fraction of UV continuum photons is 
g-TduBt^ as for a screen model, consistent with the thin 
shell approximation we make here. On the other hand, 
we predict the Lya properties of galaxies using the grid of 
models of Schaerer et al. (20111 which contains the results 
of a large number of Lya RT simulations in expanding 
shells, as described above (see [Verhamme et al.|[200^ for 
more details). The grid provides us with the Lya escape 
fraction and the emergent line profile 4?^'^“^(Lj/a) for 
more than 6000 quadruplets of parameter values (14xp, 
h, Tdust, Nhi). We interpolate the shell parameters pre¬ 
dicted by GALICS onto the grid so as to compute /esc 
and 4>(Li/a) for each individual galaxy. The observed 
Lya luminosities are thus given by I/Lya = LEya x /esc. 

A main result of this modelling is that /esc is of 
the order of unity for galaxies with a low star formation 
rate (SFR < lM©yr“^) because they are predicted to 
have low gas column densities and dust opacities (see sec¬ 
tion 3.2.3 in [Paper I| ). This suggests that the observed 
Lya luminosity could be used on average as a tracer of 
the SFR for faint LAEs. On the other hand, /esc is greatly 
dispersed between 0 and 1 for high-SFR galaxies because 
star forming galaxies can have a wide range of Nm and 
Tdust values in our model. For these galaxies, Lya is no 
longer a reliable tracer of the star formation rate. 

Our model only assumes the internal attenuation of 
the Lya line by dust in the shell, and the effect of the 
intergalactic medium is neglected. We will discuss this 
choice in more details in Section [T2| 


3 ABUNDANCES OF LBGS AND LAES 

In this section, we present the UV and Lya luminosity 
functions at z = 3-7, a redshift range that corresponds to 
the post-reionisation epoch where most dropout galaxy 
and Lyman-alpha emitter surveys have been conducted 
so far. 


3.1 UV luminosity functions 

Observationally, rest-frame UV luminosity functions (LF) 
are now rather well constrained up to z ~ 6. At higher 
redshift, significant scatter in the data remains, mainly 
due to smaller statistics, larger cosmic variance effects. 
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Figure 1. UV luminosity functions at z Ri 3, 4, 5, 6, 7 and 7.5. The solid (dashed) line is the model with (without, i.e. intrinsic) dust 


|Thompson 

l|2006 

(blue asterisks), |Gabasch et aL|||2004|l (green crosses), 

Bouwens et al.|||2007|l (red triangles), |Iwata et al.|( 

2007f 

(black asterisks), 

McLure et aL|||2009|) (blue diamonds), |Castellano et al. 

<|2010|l (light green triangles), |Bouwens et al. 

poll 

1 (red 


(|2010[ |2012[| (purple triangles). 


and more numerous contamination by foreground sources. 
In Figure we show the UV LFs at z ~ 3, 4, 5, 6, 7 
and 7.5. In each panel, the solid black line corresponds 
to the LF with dust attenuation, and the dashed curve 
gives the intrinsic LF, i.e. without dust attenuation. The 
observations, represented by the symbols with error bars, 
come from various surveys of LBGs. 


Overall, our model provides a good fit to the data 
over the redshift range of interest. While we note some 
tension at the faint end of the z = 3 LF, the agreement 
between the model and the observations is excellent at 
these magnitudes for z = 4 and 5. We somehow underpre¬ 
dict the luminosities and/or number densities of dropout 
galaxies at z = 7.5 (bottom right panel in Figure]^. Given 
that the intrinsic LF agrees with the data, our modelling 
of dust attenuation is perhaps too strong at z = 7.5. We 
do not think of this small mismatch between the model 
and the observations as a real issue since the observed 
UV LF is still uncertain at this redshift (as the fraction 
of interlopers could be quite high in z = 7-8 samples of 
dropout galaxies). 


The LBG selection probes bright UV-continuum 
sources (at ~ 1300 — 1700 A in rest-frame) and applies 
colour-colour criteria to ensure the detections at a given 
redshift, minimising the contamination by interlopers. 
The selection and filters vary from one survey to another, 
so, in our model, we choose to directly measure the 
restframe absolute UV magnitude of galaxies in an 
effective rectangular filter at 1500 ± 100 A to compute 
the UV LF (as quoted by Gabasch et al.|[2004 l, without 
any colour-colour selection. These different selections 
only introduce a little variation on the UV LFs but they 


could still be responsible for part of the slight differences 
seen in the LFs between the various surveys, and between 
the data and the model. 


3.2 Lya luminosity functions 


In comparison to UV LFs, the Lya LFs are computed 
observationally over smaller samples so cosmic variance 
effects remain significant. Moreover, most LAEs are de¬ 
tected using the narrow-band (NB) technique which in¬ 
troduces a selection in Lya equivalent width (EWQ The 
Lya LFs are therefore computed over EW-limited sam¬ 
ples that are not complete in terms of Lya luminosity, so 
caution must be taken when comparing LFs from various 
surveys which select LAEs differently (i.e. with different 
EW cuts). 

In Figure we compare the Lya LF from our model 
with available observational data at six redshifts (z = 3, 
3.7, 4.7, 5.7, 6.5 and 7). To build the LFs, we only con¬ 
sider the Lya luminosity of galaxies, and, unlike observa¬ 
tions (symbols with error bars), we do not apply any EW 
thresholds. Hence, we expect our LFs to lie above data 
points which come from LAE samples selected with high 
EW thresholds. For instance, [Ouchi et al.| ( |2008[ |2010p 
obtained four NB-selected samples of LAEs at z « 3.1, 
3.7, 5.7, and 6.6, which are often taken as a reference to 


^ All equivalent width values discussed in this paper are rest- 
frame equivalent widths. 
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Figure 2. Lya luminosity functions at z 3, 3.7, 4.7, 5.7, 6.5 and 7. The solid (dashed) line is the model with (without, i.e. 
intrinsic) dust. Symbols are observational data from [van Breukelen et al.| l |200^ (red diamonds, 2.3 < z < 4 .6), [Gronwall et al.| 
( |2007| l (brown triangles, z = 3.1) , |Ouchi et a~ l 2003 |2008| |2010[ | (green squares, z = 3.1, 3.7, 4.9, 5.7 and 6.6), [Blanc et al.| ( |201lf 
(blue asterisks, 2.8 < z < 3.8), [Rauch et ^ (20081 (orange circles, 2.67 < z < 3.75), [Cassata et al.[ p011[ | (purple downward 
triangles, z Ri 3, 4, 6), [Dawson et al.[l[2007[| (b lack crosses, z = 4.5), [Wang et ar]([2009^ (re d asterisks, z = 4.5), [Shioya et ar|([2009^ 
(black triangles, z = 4.9), [Henry et a~ |[2012[| (ora nge crosses, z = 5.7), [Shimasaku et al.[ ( [2006[ (brown circles, z = 5.7),[Hu et ab] 
( [2010[ (red triangles, z = 5.7 and 6.6), Ajiki et n)] ( [2003|[2004[[2006[ ( (light blue triangles, z = 5.7),[Murayama et al.[ ( [2007[ (orange 
triangles, z = 5.7), [Malhotra &: Rhoads[ l [2004[| (blue dot ted line, z = 6.5), [Kashikawa et al.[|[2011[| (purple dotted line, z = 6.5), 
[Shibuya et al.[([2012[l (ora nge stars, z ^ 7), [Hibon et al.| | [2012[ | (purple triangles, z ^ 7), [lye et al.[ l [2006[ l (black asterisk, z Ri 7) 
and[Vanzella et al.[|[20rT|| (green crosses, z ^ 7). 


study the statistical properties of LAEs, and notably the 
evolution of Lya LF, given the rather large number of ob¬ 
jects contained in their samples (a few hundreds at each 
redshift). In these surveys, the selection of LAEs implied 
EW thresholds which decline with increasing redshift (i.e. 
64 A, 44 A, 27 A, and 14 A at z = 3.1, 3.7, 5.7, and 6.6 
respectively), and it is unclear how this selection effect 
affects our interpretation of the redshift evolution of the 
LAE population. Especially, the EW threshold of 64 A 
at z = 3.1 coincides with the typical EW value predicted 


for galaxies with constant star formation (Chariot & Fall 


19931. Then, a significant fraction of LAEs can possibly 


be missed with such observational threshold. 

From Figure we first note that the overall agree¬ 
ment is quite acceptable given the complexity of the mech¬ 
anisms involved in the Lya transfer, and the simplicity of 
the physical picture we assume in this paper. Especially, 
our LFs compare very well to the observed ones at z « 
3.7 and 4.7 (top center and right panels in Figure]^. At 
z « 3 (top left panel), although our model agrees well 
with the spectroscop ic data of Blanc et al.| ( |2011| and 


01 


van Breukelen et al.|(2005l at Lnya > 10'‘^ erg s ^ (blue 


asterisks and red diamonds respectively), it overpredicts 
by a factor of ~ 1.5-2 the number density reported by 
[Ouchi et al.| | |2008| green squares). As already discussed 
in [Paper I[ our LF can be reconciled with their data if we 
impose high EW cuts (EWuya > 50 A provides the best 
match). 


At the faint end of the LF, our predictions seem 
to favour the number density reported by [Rauch et al.| 
(20081. These observations are however in slight disagree¬ 
ment with the findings of [Cassata et al.| ( |2011| ) who mea¬ 
sure a density of LAEs a few times lower for a similar 
detection limit (I/Lya 10"^^ erg s“^ at z « 3). Cosmic 
variance effects, due to the rather small and elongated 
volumes that are probed, as well as incompleteness issues 
or slit losses may explain the discrepancy between both 
measurements. Larger homogeneous datasets are there¬ 
fore still needed to better constrain the Lya LF at the 
faint end. This will be one of the key objectives of forth¬ 
coming instruments like the Multi Unit Spectrograph Ex¬ 
plorer (MUSE; Bacon et al.|2006 l which recently started 
to operate at VLT. We will address these issues in more 
details in a next paper (Garel et al. in prep). 

At z « 5.7 and 6.5, the EW cuts employed in NB sur¬ 
veys are small (« 15 — 25 A) so the impact on the selection 
of LAEs should be minor. As can be seen from Figure 
the shape and the normalisation of the observed Lya LF 
at z ~ 5.7 and 6.5 varies significantly from one survey 
to another, and our model agrees better with the lower 
end of the envelope of data points. For instance, at z « 
6.5, while our model reproduces nicely the observed LFs 
of |Hu et al.j (|2010| red diamonds in Figure [2|, it under¬ 
predicts by a factor of 2-4 the number densities of|Ouchi 


et al. 


(2011 


(2008 2010 green squares) and [Kashikawa et al. 


purple dotted line). The LF of Hu et al. (20101 
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was computed from their spectroscopic sample and they 
find that only « 50 % of the LAE candidates could be 
confirmed by spectroscopy. To explain the difference be¬ 
tween the LFs, |Hu et ah] (|2010[) argne that the photo¬ 
metric samples of Onchi et al. (20101 might contain high 
fractions of interlopers. However, [Kashikawa et al.| ( |201l| ) 
report that the rate of contamination in the photomet¬ 
ric sample of Taniguchi et al. (20051 at z ~ 6.5 is less 


than 20%, and their LF is in much better agreement with 
the res ults of [Onchi et al.|(|2010[ | than those o f |Hu et ah] 
(20101. Kashikawa et al. (20111 claim that the difference 


may come from the lack of completeness at the faint end 
in the sample of |Hn et al.| ( |2010p , dne to the shallow depth 
of their spectroscopic follow-up survey. 

Therefore, we conclude that our model is in agree¬ 
ment with observations at z > 6 only when comparing 
with the data that report the lowest densities of sources 
at the bright end of the LF. A weaker of effect of dust, 
or higher intrinsic Lya luminosities would be required to 


( [MT| . 


match the data of Onchi et al. (2010 \ or Kashikawa et al. 


We note that we do not take the effect of the IGM into 
account in our model whereas it is well-known that neutral 
hydrogen atoms can scatter photons on the blue side of 
the Lyo? resonance off the line of sight. This could strongly 
reduce the transmitted Lya flux, especially at z > 6. How¬ 
ever, as already shown by e.g. |Santos| (|2004[) and |Dijk-| 


stra & Wyithe (20101, Lya radiative transfer through gas 


outflows can Doppler-shift Lya photons towards longer 
wavelengths and then considerably reduce the impact of 
IGM. In Paper I (Section 4.4), we used the prescription of 


|Madau| ( 1995P to compute the mean contribution of the 
Lya forest as a function of redshift. We showed that the 
IGM had a negligible impact on the Lya luminosities of 
galaxies in our model up to z « 5 because of the peak of 
the Lya lines being redshifted sufficiently away from line 
centre by the scattering in the expanding shells. Here, we 
did a similar test, and we found that the Lya LFs up 
to redshift 7 remain unaffected by the IGM transmission, 
for the same reasons as explained above. We also tested 
the recent model of |Inoue et al.| ( |2014[ ) which predicts a 
slightly higher (lower) Lya transmission at z < 4.7 (z > 


4.7) than the recipe of Madau (19951, and gives a better 
fit to the observational data. Again, as the Lya lines of 
most our galaxies are Doppler-shifted, this IGM attenu¬ 
ation model has also almost no effect on the Lya fluxes 
at all redshifts. Nevertheless, we note that [Laursen et aL] 
(20111 have shown that the IGM can non-negligibly re¬ 


duce the flux on the red side of the Lya line at z > 6. 
Moreover, additional contributions to the IGM opacity 
might be important at z > 6, such as (i) the potential 


high number density of Lyman-Limit systems (Bolton & 
|HaehneIt]|2012[ ), or (ii) the increasing neutral fraction of 
the diffuse component of the IGM, xhi, before reionisation 


is fully complete (see Dijkstra 2014 for more details). 
Thus, a more refined IGM model with full RT treatment 
would be needed to assess the exact effect of the IGM on 
the Lya lines in our model. 

The end of the epoch of reionisation is still a sub¬ 
ject of debate, and the value of xhi at z = 6-7 is not 


fully constrained yet (e.g. Dijkstra 20141. In addition, 
Lya observations are often used as a probe of xm, e.g. 


by measuring the evolution of the Lya LF (Ouchi et al. 


20101. In this context, it is interesting to note that our 


model intrinsically predicts a decrease of 50 % of the num¬ 


ber density of bright LAEs at I/Lya = 10^® erg s ^ be¬ 
tween z = 5.7 and 6.5 without accounting for the effect of 
IGM (Figure]^. 

We now turn our interest to the cross 
UV/Lya properties of the LBG and LAE popula¬ 
tions to study their link with one another. Owing to 
different methods of selection, the dropout and narrow 
band techniques do not necessarily probe the same 
galaxies, but both populations are a subset of a common 
parent population, and it is worth asking how they 
overlap. 


4 UV/Lya PROPERTIES OF LAES 
4.1 The UV LFs of LAEs 

Here, we compute the UV luminosity functions of LAEs 
selected using Lya luminosity and equivalent width 
thresholds similar to various narrow band surveys (Fig¬ 
ure [^. In practice, we allow the observational cuts to 
vary by < 20 % in order to improve the agreement 
with the data. This is justified by the fact that the 
quoted Lya luminosity and equivalent width thresholds 
do not exactly mimic the actual colour-magnitude selec¬ 
tions ( Ouchi et al.|200'8 Dijkstra fc Wyithe|2012 I. 

Our model reproduces the UV LF of LAEs from 
[Ouchi et al.| ( |2008| ) and [Ouchi et al.j ( [2003| ) at z « 3.7 
and 4.9 respectively using selections similar to those re¬ 
ported by these authors (top centre and top right pan¬ 
els in Figure [sjl. At z > 6, our model reproduces cor¬ 
rectly the shape of the LFs observe d by [Ouchi et aL 
(20081, Shimasaku et al. (20061 and Kashikawa et al. 


(20061 but the overall normalisation is somehow too low. 


Varying the Lya luminosity and equivalent width thresh¬ 
old has little effects because the cuts are low enough such 
that few galaxies with M 1500 < —18 are removed by the 
Lya selection. As for the comparison of our Lya LFs with 
the data from the same surveys, it is plausible that the 
number density measured from photometric samples in 
these observations is overestimated due to non-negligible 
fractions of contaminants. 

At z « 3.1 however, the observed LF is strongly un¬ 
derestimated if we use the equivalent width threshold of 


Ouchi et al. (20081, as shown by the solid histogram in the 


top left panel. The value of 64 A assumed by these authors 
is close to the median Lya equivalent width in our model, 
hence a large fraction of galaxies are below this selection 
limit. We need to decrease the EW threshold from 64 A 
to 50 A to bring the model into agreement with the data 
(dotted line in the top left panel of Figure]^. In other 
terms, at a given UV magnitude, the Lya luminosities 
predicted by our model are not high enough compared to 


the LAEs of Ouchi et al. (20081. In the next section, we 


investigate the correlation between UV and Lya emission 
in more details, and we quantitatively discuss plausible 
causes for the missing high EWs in the model. 


4.2 Expected relation between intrinsic 
Lya and UV emission 

Assuming that star formation is at the origin of all the 
Lya and UV emission of galaxies, the relation between 
SFR and the intrinsic Lya lu minosity writes Z / Lva = 1.2 x 
1O'‘^(SFR/M0yr“^) erg s“^ ( Kennicu'tt|[ 1998 Osterbrock 
fc Ferland|2006 (. Following Madau et al. (19981, the rela¬ 
tion between SFR and UV luminosity density evaluated 
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Figure 3. UV luminosity functions of Lyo-selected galaxies at z 3, 3.7, 5, 6, and 6.6. The solid lines correspond to the model 
with dust attenuation included and the red filled circles are observational data l |Ouchi et al.|[^08[ |2003[ [Shimasaku et ar]|2006| 
[Kashikawa et al.|2006| l. We compared our model to each dataset individually, applying selections in terms of Lya luminosities and 
equivalent width so as to mimic real observat ions (see legends). The dashed line in the top left panel corresponds to a selection 


where the EW threshold of 
agrees quite well with the data. 


Ouchi et al. 


{2008 \ has been decreased from 64 to 50 A. With such a small change 20%), our model 


at 1500 A is Ta uv = 1-4 x 1O^°(SFR/M0yr ^) ergs ^ 
A-F 

Combining these formulae, Lhya and Misoo scale 
with one another as follow^ 

Mi5oo = -2.5 log,o ( - 18-3 (1) 

In this framework, galaxies with intrinsic luminosities 
of I/Lya = 10"^^ and 10'*® erg s“^ should display 

Mi 5 oq = -15.8,-18.3 and —20.8, corresponding to star 
formation rates of ~ 0.08,0.8 and 8MQjr~^ respectively. 

In Figure we investigate the relation between 
Lya luminosity and UV magnitude at z « 3. In the upper 
panel, the distribution of galaxies from the model, repre¬ 
sented by the black contours, shows the relation between 
Lya luminosity and UV magnitude after the effect of dust. 
Overlaid are measurements from two sets of narrow-band 
selected LAEs. The cloud of data points from [Gronwall] 
et al. ( 2007| ) (blue diamonds) covers a similar region to 
the model above their detection limit (~ 1.2 x 10"^^ erg 
s“^; blue dotted line). We note however that some data 
points lie above the model, i.e. their Lya luminosity is 
larger than predicted by the model for a given UV mag¬ 


nitude. Now comparing with data of Ouchi et al. (20081 


(red crosses), we see this time that many LAEs lie above 
the model predictions. These objects belong to the tail 
of the Lya equivalent width distribution, that is EW > 


^ These scaling relations assume a Kennicutt IMF, solar metal- 
licity, the production of hya photons through Case B recombi¬ 
nation at T= lO'^K, and are valid for constant star formation. 


100 A, which our model does not reproduce. Moreover, 
we note that most LAEs from the survey of|Ouchi et aL| 
( |2008[ ) are above the ones of [Gronwall et ah (20071 in 
the I/Lyn-Afisoo plane, which is due to very different EW 
thresholds used to select LAEs (20 A vs 64 A). 

In the middle panel of Figure we now show the 
model distribution of galaxies before applying the effect 
of dust. To show how the Lya and UV intrinsic luminosi¬ 
ties relation are expected to scale with one another as a 
function Lya equivalent width, we overplot the expected 
intrinsic Li('ya-M}|oo relation for various EWs (as labelled 
on the plot), according to: 


EWexp 

-bA,UV VXLya / 


( 2 ) 


where /3 is the UV slope used to renormalise the contin¬ 
uum luminosity density around the Lya line. Here, we 
have assumed P = —1.5, which is a typical value for 
observed LAEs (e.g. jVenemans et al.||2005| [Blanc et al.| 

[HD- 

From the bottom panel of Fig ure [4[ it is obviou s that 
many sources from the samples of Gronwall et al. (20071 
and [Ouchi et al.[ ( [2008[ ) have observed EW larger than 
in the model, even before taking the effect of dust into 
account (contours). At the same time, we note that the 
relation between intrinsic Lya luminosity and UV mag¬ 
nitude agrees well with the formula for constant star for¬ 
mation (SF) from Eq. [^(dashed green line). The produc¬ 
tion rate of Lya photons, dominated by short-lived hot 
stars, traces the very recent star formation, so the intrin¬ 
sic Lya luminosity, is proportional to the star for¬ 

mation rate averaged over the last « 10 Myr, SFR(t « 10 
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8 Garel et al. 




Figure 4. Distribution of LAEs in the Lya luminosity ver¬ 
sus UV magnitude plane at z Ri 3. The black contours show 
the number distribution of objects in the model (without any 
Lya luminosity or EW cut). Top: Lyo luminosity versus UV 
magnitude (with dust attenuation included). The data of |Gron-1 
[wall et al.| ) |2007[ l (blue diamonds) and |Ouchi et al.| i |2008[ l 
(red crosses) are also shown. The blue dotted line show the 
Lya detection limit of |Gronwall et a'r] ( |2007^ . Middle: Intrinsic 
Lya luminosity versus intrinsic UV magnitude (without dust 
attenuation included). Overlaid are the expected relations for 
various Lya equivalent widths, according to Eq.[^ Bottom: In¬ 
trinsic Lyo luminosity versus intrinsic UV magnitude (without 
dust attenuation included). The green dashed lines show the 
expected relation in the case of constant star formation (Eq. 

and in the case of an instantaneous starburst (see text). 
Same data as above (i.e. not corrected for dust). 


Myr) (Chariot & Fall 19931. The UV continuum intensity 
is however a direct function of the mean SFR over longer 
timescales, of the order of a few hundreds Myr (La,uv oc 
SFR(t ~ 100 Myr); [Madau et al.||199^ . In the constant 
SF scenario, i.e. SFR(t « 10 Myr) = SFR(t rs 100 Myr), 
the Lya equivalent width reaches an equilibrium value of 
approximatively 60 A after « 10 Myr, which is consistent 
with the median intrinsic EW value we find in our model. 

In GALICS, star formation is iteratively computed at 
each timestep using the Kennicutt law, namely Esfr oc 


and a gas surface density threshold: Eg 


= 10 ^ 


observations, is in broad agreement with the Toomre 
model for gravitational instability-triggered star forma¬ 
tion ( |Kennicutt||1989| . In GALICS, this criterion is met 
most of the time in LAEs and LBGs we consider in this 
paper, hence the star formation appears to be close to a 
constant process, although there is scatter in the 
MJsoo relation due to small variations of the recent star 
formation history or stellar metallicity. We al so consider 
here the LiJ'ya-AfJsQo relation (see Fig. 15 of Verhamme 
|et al.|[^008[ ) that is expected a few Myr after a burst of 
star formation (dashed green line). We see from the bot¬ 
tom panel of Figure |4|that the high-EW sources of|Ouchi| 
[eTaL] ( (20081 ) can be interpreted by this scenario, which 
produces higher Lya luminosity for a given UV emission 
than ongoing constant SF. Thus, the high-EW LAE pop¬ 
ulation probed by [Ouchi et al.| ( |2008[ ) could correspond 
to starburst galaxies, rather than constantly star forming 
objects. 

Although our model can recover the standard 
Lya equivalent widths (EW ~ 0-70 A) which seem to 
make the bulk of the LAE population ((Cassata et al.| 
[MT| , it certainly requires extra ingredients to explain 
the high-EW galaxie^ 


4.3 High Lya equivalent widths 


As discussed in the previous section, the commonly ob¬ 
served of high-EW LAEs (EW > 100 A) are missing in 
our model, even when considering intrinsic values. Here, 
we discuss various potential causes responsible for this 
mismatch. 

First, it is often claimed that Lya radiation trans¬ 
fer effects could enhance the intrinsic EW, as predicted 
by the analytic model of Neufeld ( 1991[) for a mult iphase 
medium. (Laursen et al.| ( |2013[ ) and Duval et al.| ( |2014[ ) 
have shown that the (angle-averaged) boost of Lya EWs 
in a clumpy multiphase medium would require very spe¬ 
cial conditions unlikely to exist in most high-redshift 
galaxies. However, the recent model of |Gronke fc Dijk-| 
(2014) suggests that significant Lya enhancements 


stra 


can be found along particular sight lines where UV con¬ 
tinuum photons are absorbed by dust clouds, whereas 
Lya photons escape isotropically. Alternatively, [Vct^ 
(hamme et al. ] ( |M^ have reported a strong inclina¬ 
tion effect using Monte-Carlo Lya RT in hydrodynami- 
cal high-resolution simulations of galactic disks, in which 
the Lya equivalent width measured perpendicularly to 
the disk is larger than in the edge-on direction. While 
both UV continuum and Lya are emitted isotropically in 
the disk, a fraction of the Lya photons emitted edge-on 


^ Here, we have focussed on z = 3 LAEs, but the issue is similar 
at higher redshifts, as large EW values are also measured in 
LAEs up to z Ri 6 (e.g. [Kashikawa et al.|2011||. 
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The UV, Lyman-a, and dark matter halo properties of high redshift galaxies 9 


may be scattered off by hydrogen in the face-on direc¬ 
tion, therefore boosting the observed EW along the face- 
on line-of-sight. 

Second, high EWs can be intrinsically produced in 


galaxies. For instance, Schaerer (20031 show that the EW 


strongly depend on the metallicity of stars and the IMF. 
Our SED libraries ( [Devriendt et al.]|1999| ) only take into 
account metallicities larger than Z = 2.10“^ so we cannot 
investigate the impact of metal-poor stars in our model, 
but we can compare the EW variation between the re¬ 
sults from our fiducial Kennicutt IMF with a low-mass 
cut-off of 0.1 Mq (labelled IMFa) and from a more top- 
heavy Kennicutt IMF that assumes a low-mass cut-off of 
4 M 0 (IMFb). In Figure]^ (top panel), we show the z = 3 
model distributions of intrinsic EWs (before dust attenu¬ 
ation included; solid histogram), and observed EW (after 
dust attenuation included; dashed histogram) with IMFa. 
This model is in reasonable agreement with the data of 
|Gronwall et al.| ( |2007[ ) (blue crosses) in the range 20-70 
A, which corresponds to the peak of the EW distribu¬ 
tion. However, the high-EW tail extending to ~ 250 A is 
poorly matched, as expected from the above discussior]^ 
IMFb tends to produce higher EWs due to a larger ratio 
of high-mass to low-mass stars (orange dotted curve), but 
the distribution remains as narrow as for IMFa, and it 
does not reproduce the high-EW tail either. Alternative 
IMFs could have a stronger impact on the EW distribu¬ 
tion, as shown by |Orsi et ah] ( |2012[ ) who use an extreme 
top-heavy, fiat IMF. They find a much broader EW dis¬ 
tribution than ours, but it turns out to be even too broad 
compared to observations, and it peaks to higher EW than 
observed (see their Figure 11). We also note that, rather 
than invoking very top-heavy IMFs, [Forero- Romero &: Di-| 
|jkstra| | |2013[ ) demonstrate that the stochastic sampling of 
the IMF can lead to strong variations in terms of EW, 
and significantly broaden the EW distribution. 

As mentioned in Sec. |4.2| a starburst should be able 
to produce higher EWs than ongoing star formation at a 
constant rate, as it is the case in GALICS. We test this 
scenario in the bottom panel of Figure where we com¬ 
pare the intrinsic EW distribution from our hducial model 
(black solid histogram) with a modified version of GAL¬ 
ICS in which star formation has a stochastic duty cycle of 
10% (’stochlO’; red dashed histogram). In GALICS, star 
formation (and all other galaxy formation processes) are 
integrated over timesteps of 1 Myr. Here, we randomly 
switch on star formation once every ten timesteps to see 
the effect of bursty star formation onto the Lya equivalent 
widths. With this ad-hoc model, we clearly see that the 
EW distribution is broader than the fiducial model’s dis¬ 
tribution, and agrees much better with the data of |Gron-] 
[wall et al.| ( |2007| ). 

An obvious method to mimic a SF duty-cycle is to 
set a high gas surface density threshold. While the cri¬ 
terion = 10 ^^ cm“^ is nearly always met in our 

fiducial model for galaxies considered here, we tried to in¬ 
crease this value to 10^^ cm“^. In this case, accreted gas 
needs to accumulate for longer periods in the galaxy in 
order to reach the surface density threshold, at least in 
some objects. The corresponding intrinsic EW distribu¬ 
tion plotted as a green dot-dashed histogram on the bot- 


^ It is worth pointing out here that our simulation box is few 
times larger than current LAE surveys, so it is unlikely that 
these very strong emitters are missing in our model because of 
our finite volume. 
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Figure 5. Lya equivalent widths distribution of LAEs at z 
= 3. Top: The black solid and dashed histograms show the 
EW distributions of LAEs predicted by our model before and 
after the effect of dust included, respectively. Our fiducial 
model assumes a Kennicutt IMF, with a low-mass cut-off of 
0.1 M 0 (IMFa). For comparison, we also show the EW distri¬ 
bution (without dust) using a more top-heavy Kennicutt IMF, 
with a low-mass cut-off of 4 M 0 (IMFb; orange dotted curve). 
LAEs have been selected using Li^ya ^ 1-2 x 10*^^ erg s~ 


Gronwall et al. 


(2007 1 (blue 


EW ^ 20 A, as in the survey of 
crosses). Bottom: To highlight the effect of starbursts on the 
Lyo equivalent widths, we compare the intrinsic EW distri¬ 
bution from our fiducial model (solid black line), with similar 
models in which (i) star formation occurs stochastically with 
a duty cycle of 10 % (dashed red line), and (ii) star formation 
only occurs when the gas surface density is larger than 10 ^^ 
cm“^ in galaxies (dot-dashed green line). 


tom panel of Figurej^is similar to what we find for our ad- 
hoc ’stochlO’ model, with an extended tail towards high 
large EW valnes. Although star formation occurs more 
rarely in these two models than in our fiducial one, we 
note that the intrinsic Lya and UV LFs (and the stellar 
mass functions) are not changed by a large amount. In¬ 
deed, more gas accumulates in the galaxy for the ’stochlO’ 
and high-Egag'*®*' models, so SF events turn more gas into 
stars at once using the Kennicutt law. These models even 
produce slightly larger SFRs, hence larger Lya and UV 
luminosities, due to the power-law index of the Kenni¬ 
cutt law being greater than 1. However, the higher gas 
content would increase the effect of dust attenuation, as 


both quantities are correlated (Section 2.2 in Paper U, so 
these models will no longer necessarily give a good match 
to the observed UV and Lya LFs. Overall, stochastic star 
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formation scenarios constitute a viable mechanism to pro¬ 
duce high Lya EWs in high redshift galaxies, and we will 
investigate these models in more details in future studies. 


5 UV/Lya PROPERTIES OF LEGS 
5.1 The fraction of emitters within LBGs 


Spectroscopic follow-ups of LBG samples show that 
Lya emission is more often detected in UV-fainter galax¬ 
ies. In Figure]^ we test our (fiducial) model against the 
observations of [Stark et al. ( |2010| who found that the 
fraction of strong Lya emitters (EW > 50 A), -Ylae, in¬ 
creases from Ri 7 % (20 %) at Misoo = —22 to 45 % (55 
%) at Misoo = —19 for z « 4 (z ~ 5). We find a very 
similar trend at both z « 4 and 5.5, as shown by the 
histograms in the top and bottom panels. The agreement 
with the data of Stark et al. (20101 (represented as circles) 
is excellent at Misoo < —19.5 for z ~ 4. At z ~ 5.5, the 
fraction of LAEs with EW > 50 A in the model is slightly 
lower than the observations. However, as shown by the 
red dotted histogram, we find a good match to the data 
if we decrease the EW threshold by only 10 % (EW > 45 




Figure 6. Fraction of strong Lyo emitters among samples of 
LBGs at z Ri 4 (top) and 5.5 (bottom). In each panel, the 
histogram shows the predicted fraction of LAEs Alae in five 
bins of rest-frame absolute UV magnitude Misoo. The model 
is compared to the data of jStark et al.| ( |2010^ (represented by 
circles with error bars) who measured the fraction of LBGs 
with Lyo equivalent widths larger than 50 A. As shown by 
the dotted line in the bottom panel, we find a slightly better 
agreement with the observations at z Rs 5 with an EW cut of 
45 A instead of 50 A (solid line). 



20.0 20.5 21.0 21.5 22.0 22.5 

log Nhi [cm-'] 


Figure 7. Lyo/UV escape fraction ratio in z = 3 LAEs as a 
function of Hi shell column density. The dots represent indi¬ 
vidual galaxies in our model, colour-coded relatively to their 
Lya EW (after applying dust attenuation). Eor comparison, we 
show the Lya/UV escape fraction ratio expected from constant 
Lyo escape fraction models | |Le Delliou et al.|2005||Nagamin^ 
jet al.|2010|l, assuming fesc(Lyo) = 20% (contours). 


A). The fractions predicted by our model in the faintest 
bin are a bit lower than their observational measurement, 
but the error bars are significant at this magnitude. 

The Xlae evolution as a function of Misoo is related 
to the anti-correlation between Lya equivalent width and 
UV magnitude which is commonly observed in LBGs and 


LAEs ( Ando et al.|2006 Stanway et al.|2007 Ouchi et al. 
|2008| ). As discussed in [Paper I| we find that UV faint 
galaxies are slightly more likely to display high intrin¬ 
sic EW than bright LBGs. However, the main driver of 
the observed Alae-Misoo relation in our model remains 
the differential dust attenuation experienced by Lya and 
UV photons, caused by Lya resonant scattering in Hi 
gas. This is illustrated on Figure where we plot the 
ratio of UV to Lya escape fractions as a function of 
the shell Hi column density, Nhi, for galaxies at z = 3 
which have I/Lya > 10'*^ erg s“^. Each dot represents a 
galaxy, colour-coded as a function of its Lya equivalent 
width. The fesc(Lya)/fesc(UV) ratio is approximatively 
one at low column density, but it sharply decreases to¬ 
wards larger values, which highlights the stronger effect 
of dust on Lya due to resonant scattering in Hi optically 
thick media. High-SFR galaxies (i.e. UV-bright galaxies) 
having a larger Hi column density on average in GALICS, 
the Nui-dependent feac(Lya)/fesc(UV) ratio unavoidably 
leads to a lower Lya fraction in bright LBGs compared 
to faint LBGs (Fig. |^. 

For comparison, we have overlaid on Figure the 
predictions for the constant Lya escape fraction model 
(contours), which is often used in the literature to fit 
the observed Lya luminosity functions (e.g. [Le Delliou[ 
et al.|200^ Nagamine et al.|2010 1 . Fixing fesc(Lya)=20%, 


which is the value needed to roughly match the bright-end 
of the z = 3 LF in our model, the feac(Lya)/fe 3 c(UV) ratio 
increases towards large Nm values. A consequence of this 
trend, opposite to what is found with our fiducial model, 
would be to predict a higher LAE fraction in brighter 
LBGs, at odds with the observations. 

The scaling of the ratio of UV continuum escape frac¬ 
tion to resonant Lya escape fraction with Nhi and SFR, 
is then a key factor to explain the Alae-Misoo trend, as 
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also reported by |Forero- Romero et al.| | |2012[ ), who coupled 
hydrodynamical simulations with Lya RT in galaxies ap¬ 
proximated as dusty slabs. 


5.2 The evolution of the Lya fraction in LBGs 
with redshift 


According to recent studies, the fraction of strong 
Lya emitters among LBGs, Xlae, seems to increase from 
z « 3 to 6 ( [Stark et al.pOlTj jCassata et al.pOlS l. This 
trend appears to invert at z = 7, as reported by jOno et al.j 
(20121 and Schenker et al. (20121 who find that Xlae dra¬ 


matically drops between z « 6 and 7 (see also jPentericcij 
et al.|2011 l. The redshift evolution of Xlae could be due 
to a rise of the mean internal Lya escape fraction from 
galaxies at z « 3-6, and by a sharp increase of neutral 
IGM opacity at z > 6. Nonetheless, it remains plausi¬ 
ble that the apparent trend of the fraction of emitters in 
LBGs with redshift is also (i) affected, or even driven, by 
the intrinsic evolution of the Lya, UV, ionising properties 
of galaxies ( Dayal &: Ferrara]|2012 Dijkstra et al.|[2014 |, 
(ii) due to a reduction of the Lya transmission caused by 
the increase of the number of Lyman Limit systems at 
the end of the epoch of reionisation ( Bolton fc Haehneltj 
20121 , and/or (iii) altered by small number statistics (see 
Dijkstra]|2014| [Mesinger et al.||2015[ for a more detailed 
discussion). 

Even though the study of reionisation is out of the 
scope of this paper, it is interesting to compare our predic¬ 
tions to the observed redshift evolution of Xlae without 
taking the effect of IGM into account. Figurej^shows the 
evolution of Xlae of faint (in red) and bright LBGs (in 
black). The symbols are observational estimates of Xlae 
where a LAE is defined as having EW > 55 A in the top 
panels, and EW > 25 A in the bottom panels. The solid, 
dashed, dotted and dot-dash curves are predictions from 
our model for several EW cuts: EW ^ 25 A, EW > 35 
A, EW ^ 45 A, and EW ^ 55 A respectively. In all pan¬ 
els, we clearly see that the LAE fraction is reduced at 
all redshifts when higher EW cuts are applied. There is 
a decreasing trend with redshift for EW ^ 25 A, but it 
flattens for higher EW cuts. 

Erom the upper panels (EW > 55A), we see that we 
need to lower the EW cut from 55 A to 45 A to match 
the data (dotted lines). In this case, the Alae computed 
from the model is based on small statistics because there 
are few LBGs in the corresponding ranges of magnitude 
with EW ^ 55 A, so it is hard to draw robust conclusions 
from our model. For the EW ^ 25 A measurements how¬ 
ever (bottom panels), a better agreement is obtained if 
we increase the threshold to select only emitters stronger 
than EW = 35 A. 

In order to clarify the disagreement between the 
model and the observed redshift evolution of Alae in 
LBGs, as well as the strong variation with the EW cuts, 
we now investigate the variation of the Lya equivalent 
width distributions with UV magnitude. Indeed, by con¬ 
struction, should the model perfectly match the observed 
EW distributions of UV-selected galaxies, the redshift 
evolution of the LAE fraction in LBGs would be recov¬ 
ered. These are shown in Figure]^ for z = 3, z = 4, and z 
= 5.5 as a function of Misoo, and compared to the VUDS 
data of [Gassata et al.| ( |2015| ) (green crosses). 

First, we note that the model distributions shift to 
low EW values as brighter and brighter LBGs are consid¬ 
ered. They span a range similar to the observed ones, from 



EW>55A 


EW>25A 


34567893456789 
z z 


Figure 8. Fraction of Lya emitters, Xlae? a function of 
redshift among samples of faint LBGs (—20.25 < M 1500 < 
— 18.75; left panels, in red) and bright LBGs (—21.75 < 
^1500 < —20.25; right panels, in black). We show the observed 
fraction of strong Lyo emitters (EW ^ 55 A) and weaker emit¬ 
ters (EW ^ 25 A) in the top and bottom panels respectively. 
Data: Crosses ([Gassata et al.||2015| l, dia monds ([Stark et al.| 
|2011[[Schenker et al.|2012| |, and triangles ( jOno et al.|2012| l. In 
all panels, we compare the observations to the LAE fraction 
predicted by our model for various EW cuts (solid line: EW 
^ 25 A; dashed line: EW ^ 35 A; dotted line: EW ^ 45 A; 
dot-dashed line: EW ^ 55 A). 


« -40 A to 70 A. Unlike LAE surveys, which preferentially 
pick up low-continuum sources, UV-selected galaxies show 
very few high EWs (> 70 A). Interestingly, the tail of the 
distributions appears to be slightly shifted to larger EW 
values for fainter LBGs. This explains why the fraction 
of galaxies with EW > 50 A is found to increase towards 
faint UV magnitude in Eigure in agreement with the 
data of Stark et al. (20101. However, the fraction of galax¬ 
ies drops sharply at EW « 50 A for all UV magnitude 
selections considered here, echoing the rapid variation of 
Alae for Lya equivalent width thresholds near 50 A, as 
discussed in Figure 


Second, the VUDS EW distributions computed for 
galaxies with Misoo ^ -21 peak to lower EWs than 
the model distributions for the same UV magnitude 
range (green dashed curves), which is the reason why 
we find larger Alae than what is observed for EW ^ 
25 A (bottom panels of Fig ure [s i. However, the EW 
distributions of Gassata et al. (2015 I agree very well with 
the brightest LBGs in the model (Misoo < -22.5), as 
shown by the red dot-dashed curves. 


Although our model catches well the trend between 
Lya equivalent widths and UV magnitude, the above 
study suggests that the exact scaling between these quan¬ 
tities is very sensitive to the cuts used to compute the 
Lya fractions in LBGs. In addition, some fine-tuning 
of the model is required to fully reproduce the ob¬ 
served UV/Lya cross-properties. In particular, the ratio 
of UV/Lya escape fraction is found to have the right scal¬ 
ing with relevant physical properties, like Hi column den¬ 
sity, or SFR, in order to interpret the link between the 
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EW [A] EW [A] EW [A] 

Figure 9. Lya equivalent widths distributions of UV-selected galaxies at z = 3 (left), z = 4 (center), and z = 5.5 (right). In each 
panel, we show the distributions for different UV magnitude cuts, as labelled. The green crosses correspond to the data of |Cassata| 
|et al.| < |2015[ l (rebinned) for galaxies selected with M1500 5 -21 (see their Figure 1). 
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Table 1. Properties of LAEs and LBGs at z = 3, 3.7, 5, 6, 6.6 and 7.5 for three bins of Lyo and UV luminosity. Here, LAEs are 
defined to have EW > 0 A and they are classified as faint (10^^ < L-Lya < 10^^ erg s 1), typical < Li^ya < 10'*^ erg s l), 

and bright (10^^ < Li^ya < 10"*^ erg s“^) according to their Lya luminosity. The faint, typical, and bright samples of LBGs are 
selected from their UV magnitude, —15.8 > M1500 > —18.3, —18.3 > M1500 > —20.8, and —20.8 > M1500 > —23.3 respectively. 

median host halo masses of LAEs and LBGs in each subsample. The subscripts and superscripts next to 
each value correspond to the 10th and 90th percentiles respectively. For each subsample (faint, typical and bright), the third row 
shows the ratio of the number of LAEs to the number of LBGs, n'^^® / 


observed properties of LAEs and LBGs, but a stronger 
differential LfV/Lya dust extinction in UV-bright objects 
is needed to improve the quantitative agreement with the 
data. 


6 HOST HALOS OF LAES/LBGS 

In this section, we analyse the properties of the host ha¬ 
los of LAEs and LBGs as predicted by our model. While 
the study of the clustering of LAEs and LBGs will be 
addressed in a companion paper (Garel et ah, in prep), 
here, we focus on the halo masses and halo occupation of 
LAEs and LBGs as a function of redshift to identify how 
these galaxies respectively trace the dark matter struc¬ 
tures. We also briefly discuss how the dynamical range 
probed by LAEs depends on the EW selection. 


6.1 Halo masses as a function of Lya /UV 
luminosities 

On the one hand, clustering analysis reveal that LAEs 
reside in rather low mass halos, with a median mass of 


~ at z « 3 (Gawiser et al. 20071. Overall, the 

measured biases of LAEs at z « 3-7 correspond to dark 


matter halos in the range Mh = 10^° (Ouchi et al. 


|2010[ ). On the other hand, [Hildebrandt et al.| ( |2009 1 find 
the typical halo mass of LBGs at z = 3-5 to be > 10 ^^Mq, 
namely about one order of magnitude more massive than 
those of LAEs ( Hamana et al.||2004[ [Ouchi et al.|[2004| 
|2005||McLure et al.||2009 K Nevertheless, large uncertain¬ 
ties remain in the determination of the halo masses, no¬ 
tably for LAEs for which the selection and the detection 
limit can be quite different from one survey to another. 
Moreover, samples of narrow-band selected galaxies are 
likely to contain significant fractions of interlopers, which 
lead to an underestimation of the halo mass. Correcting 
for such effect, Kovac et al. (20071 derive a similar bias 
for LAEs and LBGs at z « 4. However, they report that 
LAEs are 2-16 times rarer than LBGs for a given halo 
mass. We propose to study these considerations with our 
model in the next paragraphs. 


In Table we present the median halo masses 
in three bins of observed (i.e. dust-attenuated) 
Lya luminosity and UV magnitude at redshifts z = 
3 to z = 7.5 in our model. The Lya luminosity ranges 


© 2014 RAS, MNRAS 000,[iHT7| 







































The UV; Lyman-a, and dark matter halo properties of high redshift galaxies 13 




log M, [M,] log M, [M,] 

Figure 10. Halo occupation distribution of LBGs. The mean number of galaxies per dark matter halo, (A^g}haio’ ^ function of 
halo mass is shown by the black solid line for z = 3 (left panel) and z = 6 (right panel). The coloured curves represent the mean 
halo occupation of faint (—15.8 > M1500 > —18.3: blue dotted line), typical (—18.3 > M1500 > —20.8: green dashed line), and 
bright (—20.8 > M1500 > —23.3: red dot-dashed line) LBGs. 




Figure 11. Halo occupation distribution of LAEs. The mean number of galaxies per dark matter halo, (A^g}haioi ^ function of 
halo mass is shown by the black solid thick line for z = 3 (left panel) and z = 6 (right panel). The coloured thick curves represent 
the mean halo occupation of faint (10“^^ < Lj^ya < 10"^^ erg s“^: blue dotted line), typical (10“^^ < I/Lya < 10^^ erg s~^: green 
dashed line), and bright (10“^^ < Li^ya < 10^"^ erg s~^: red dot-dashed line) LAEs, assuming no equivalent width selection. The 
thin curves are similar but for LAEs with EW > 40 A. 


are 10"'^ < I/Lya < 10"'^ erg s 10"^^ < L-Lya < 10"^® erg 

and < L-hya < 10"^^ erg s and we term them 
faint, typical and bright LAEs. Similarly, the notation of 
faint, typical and bright LBGs will be used for galaxies 
with —15.8 > Mi5oo > —18.3, —18.3 > M 1500 > —20.8, 
and —20.8 > M 1500 > —23.3. The limit values correspond 
to SFRs of ~ 0.08,0.8,8 and 80 Mq yr“^ according to 
the expected intrinsic scaling between Z/Lya, Afisoo and 
SFR (see Section 4.2 1 . 

First, from inspection of Table we see that the 
host halo mass of LAEs and LBGs increases from faint 
to bright luminosity/magnitude at all redshifts. This is 
expected because the gas accretion rate is higher in more 
massive halos, so galaxies have higher SFR, hence higher 
intrinsic Lya and UV emission. Second, for the three lu¬ 
minosity/magnitude bins we defined, LAEs and LBGs in¬ 
habit very similar halos, in agreement with the results 
of |Kovac et al.| ( |2007[ ). For instance, at z = 3, the me¬ 
dian halo mass is S.IO^^Mq for faint LAEs/LBGs, and 
1 — 2.1O'^^M0 for bright LAEs/LBGs. This result is sim¬ 
ply due to Lya and UV emission being both intrinsically 
powered by SFR, as discussed in Sec. |4.2| 

Although we find very similar host halo masses for 
LAEs and LBGs at all luminosities considered here, their 


occupation rate of dark matter halos significantly evolves 
(see Table [^. The ratio of the number of LAEs to the 
number of LBGs, n*"^® / is about one for faint ob¬ 

jects because, at the faint end, dust extinction has little 
effect so nearly all LBGs are also LAEs, and their ob¬ 
served UV magnitudes and Lya luminosities follows the 
expected relation given by Eq.[^ At the bright end, radia¬ 
tive transfer effects affect more strongly the Lya intensity 
than the UV so that only 20 — 30% of bright LBGs will 
have a bright Lya counterpart. In Tablewe considered 
LAEs to have EW > 0 A. If we assumed a higher EW 
threshold (e.g. 10-30 A), the ratio of the number of LAEs 
to LBGs would still be about one for faint sources but 
it would be less than 20 — 30% for bright galaxies, since 
UV bright sources have lower EWs than faint LBGs on 
average (as discussed in Figure]^. 

Several reasons can be invoked to explain why LAEs 
are observationally measured to be located in less mas¬ 
sive halos than LBGs: (i) the LAEs are selected in deeper 
surveys than LBGs, (ii) the LAE samples were highly con¬ 
taminated, or (iii) the LAEs were selected in the tail of 
the EW distribution where Lya traces a very recent star- 
burst, so their Lya luminosity is higher than expected for 
a given UV magnitude (see Eq.|^. 
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Figure 12. Detectability of LAEs as function of host halo mass for various EW thresholds for z ~ 3 (left panel) and z ~ 6 (right 
panel). The ’LAE detection rate’ is the ratio of the number of LAEs with EW > 20, 35, and 50 A (dotted, dashed and dot-dashed 
curves respectively) to the number of galaxies with a Lyo emission (i.e. EW > 0 A). Here, we have only considered galaxies brighter 
than Li^ya > 10'*^ erg s“^, which corresponds approximatively to the sensitivity limit of current photometric surveys of LAEs. 
Galaxies with such luminosities only inhabit halos that are massive enough, hence the apparent cut-off of the curves at low Mj,. 


6.2 Halo occupation of LAEs and LBGs 

In order to investigate how LBGs and LAEs respectively 
populate dark matter halos, we now turn our interest to 
the halo occupation of UV- and Lyo-selected galaxies. 
Figures [To| and [TT] plot the mean number of galaxies per 
halo, as a function of halo mass, for z ~ 3 (left 

panels) and z ~ 6 (right panels). The black solid curves 
in all panels correspond to all galaxies. While our model 
predicts about one galaxy per halo at low Mh, 
increases with halo mass and reaches about 1,000 ( 100 ) 
in the most massive structures at z « 3 (z « 6 ). 

Here, we again split galaxies according to their ob¬ 
served UV and Lya luminosities: the thick dotted blue, 
dashed green and dot-dashed red curves refer to faint, 
typical, and bright objects respectively. Although galax¬ 
ies can form in halos with masses > 2 . 10 ® Mq in our sim¬ 
ulation, we can see that there is a minimum allowed host 
halo mass for a given UV or Lya luminosity threshold. For 
instance, almost no faint LBG (—15.8 > M 1500 > —18.3) 
inhabit halos less massive than 10^® Mq at z = 3 (blue 
dotted curve in left panel of Figure 101. Moreover, the 


minimum halo mass increases towards brighter objects. 
Again, this is the direct consequence of the UV luminos¬ 
ity varying as the SFR, and the SFR being proportional 
to the gas accretion rate, and thus to the halo mass. 

Each halo mass bin is populated by one correspond¬ 
ing (central) galaxy of given luminosity and by satellites 
which belong to a fainter population: there is one cen¬ 
tral faint LBG per halo of log(Mh/M 0 ) « 10.5, one cen¬ 
tral typical LBG per halo of log(Mh/M 0 ) « 11.5, and 
one bright LBG per halo of log(Mh/M 0 ) « 12.5. A sim¬ 
ilar behaviour is seen for LAEs (Eigure [TT]), but we note 
that the number of bright LAEs in massive halos is less 
than the number of bright LBGs. This echoes the results 
from Tablej^where we found that bright LAEs were rarer 
than bright LBGs (see Section |6.1| ), and this is due to 
the stronger effect of dust attenuation for Lya than UV 
continuum which can turn intrinsically bright LAEs into 
much fainter objects (see Section [4.2[ |. 

Still in Eigure [m the thin curves correspond to 
strong emitters only (EW> 40 A). For faint LAEs, 
this criterion does not remove many objects. However, 
(■^g)haio starts decreasing more significantly for the typ¬ 


ical and bright samples, especially for galaxies residing 
in massive halos. This indicates that massive halos are 
more likely to host weak emitters (low EW galaxies) than 
low-mass halos in our model. 


6.3 The bias on LAE host halos introduced by 
EW cuts 

As noticed in Figure EH selecting LAEs above a given 
Lya EW especially removes galaxies in more massive ha¬ 
los. Here, we investigate this aspect more quantitatively. 
In Figure |12[ we plot the LAE detection rate as a func¬ 
tion of host halo mass at z « 3 (left panel) and z « 6 
(right panel), adopting various EW cuts. We consider only 
galaxies brighter than 10 '*® erg s“*, and we define the 
LAE detection rate as the ratio of the number of LAEs 
with EW > 20, 35, and 50 A (dotted, dashed and dot- 
dashed curves respectively) and the number of all LAEs 
(i.e. galaxies with Lya in emission, EW > 0 A). 

We find that the LAE detection rate decreases at all 
halo masses as the EW threshold increases. Using a value 
of 20 A still allows to detect the vast majority of LAEs. 
However, many line emitters may be missed when larger 
values are adopted according to our model. Moreover, for 
a given EW cut, the drop of the LAE detection rate is 
more significant in massive halos than low-mass ones at 
both z ~ 3 and z « 6 . As an example, for EW = 35 A at 
z « 3, it decreases from ~ 90% at Mh = 1 O**M 0 to 60 % 
at Mh = 1 O*®M 0 . This behaviour has two main causes. 

First, central, intrinsically Lya-bright, galaxies which 
are located in massive halos are often strongly affected by 
dust extinction which can significantly reduce, not only 
their Lya luminosity, but also their equivalent width. Sec¬ 
ond, the bulk of galaxies residing in massive halos have a 
lower intrinsic Lya equivalent width than sources located 
in low-mass halos in our model. They are often satellite 
galaxies, so they no longer accrete fresh gas from the IGM. 
Thus, they have less intense recent star formation, which 
leads to lower intrinsic ratio of ionising (i.e. Lya) to UV- 
continuum photons. 

Therefore, the host halo population probed by 


narrow-band surveys using high EW selections (e.g. Ouchi 
|et al.|2008|) may be biased towards low-mass halos. Never 
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theless, we remind our model does not reproduce very well 
the EW distributions from narrow-band LAE surveys. As 
discussed in Section [4.2[ we underpredict the number of 
high EW sources, i.e. EW > 100-150 A, so the curves for 
the LAE detection rates in Figure should be shifted 
up by an amount which depends on the exact form of 
the EW distribution. On the other hand, these high-EW 
sources correspond to UV-faint galaxies according to ob¬ 
servations so they will be located in low mass halos on 
average. Should they be added to our model, this could 
only have the effect of accentuating the trend between 
LAE detection rate and halo mass that we see in Figure 

m 


7 DISCUSSION - SUMMARY 


In this paper, we have investigated the connection be¬ 
tween Lyman-alpha Emitters and Lyman-Break galax¬ 
ies from z « 3 to z « 7 with a semi-analytic model 
of galaxy formation. Observationally, LAEs are detected 
from their strong nebular Lya emission line whereas LBG 
surveys pick up strong near-UV stellar continuum. While 
both channels intrinsically trace the stellar formation in 
galaxies, it remains unclear how different are the galaxies 
probed by these two techniques. 

We have used the GALICS hybrid model to de¬ 
scribe the formation and evolution of galaxies, based on 
a new high resolution N-body dark matter simulation 
run with a set of cosmological parameters consistent with 
the WMAP-5 data release. The simulation box contains 
more than a billion particles in a representative volume 
of the Universe (100 h~^ cMpc on a side) which allows us 
to study the statistical properties of currently observed 
LAEs and LBGs at high-redshift. We describe the radia¬ 
tive transfer of Lya and UV photons through dusty gas 
outflows using a simple expanding shell model. To account 
for the resonant scattering of Lya photons, our model is 


coupled to the library of Schaerer et al. (20111 which pro¬ 


vides us with the Lya escape fraction and the emergent 
line profile based on 3D numerical Monte-Carlo radiation 
transfer simulations. 

With this approach, we aimed at interpreting the 
large sets of observational constraints on LAEs and LBGs 
to better understand the link between these two galaxy 
populations. We first focused on the comparison of our 
model with Lya and UV observations of Lya- and UV- 
selected galaxies. In a second part, we presented predic¬ 
tions on the host halo masses and halo occupation of LAEs 
and LBGs based on our model. 

We summarise the main results from our paper as 
follows. 


thresholds. This highlights the fact that strong emitters 
(EW > 70A) are almost absent from our model, and we 
then investigate in more details the link between Lya and 
UV emission properties of LAEs, i.e. the Lya equivalent 
widths. 

(iii) While various mechanisms, extensively discussed 
in the literature, can produce high EWs, we suggest that 
bursty star formation, rather than almost constant star 
formation, should bring our model in agreement with 
the observed EW distributions of LAEs. We show that 
increasing the surface density SF threshold by an order 
of magnitude compared to the standard threshold can 
mimic stochastic star formation rate, needed to match 
the EW distributions. 


(iv) We find that the fraction of strong line emitters 
(EW > 50 A) in LBG samples increases towards faint 
UV magnitudes, as reported by Stark et al. (20101 at z = 
4-6. This trend is essentially due to resonant Lya photons 
being more affected by dust than UV continuum photons 
in bright LBGs because these objects have a higher 
Hi column density on average, echoing the results of 


Verhamme et al. (20081. 


(v) We study the redshift evolution of the frac¬ 
tion of strong (EW > 55 A) and weak (EW > 25 A) 
Lya line emitters, Xlae, within UV-selected samples. 
Overall, Xlae seems to quickly evolve when varying UV 
magnitude limits and Lya equivalent width cuts. The 
predictions agree reasonably well with the data for strong 
emitters, but the low number statistics and the rapid 
variation of the shape/normalisation of the XcAE-redshift 
relation with UV magnitudes limits and Lya equivalent 
width cuts prevent us from drawing robust conclusions. 
In addition, as shown by the study of the model EW 
distributions as a function of UV magnitude, a stronger 
differential UV/Lya dust extinction in UV-bright objects 
is needed to improve the quantitative agreement with the 
data. 


(vi) We find that LAEs and LBGs in each sample 
are located in very similar halos, and that brighter 
sources inhabit more massive halos average. At z = 3 for 
instance, faint LAEs and LBGs are hosted by 5 x 10^° Mq 
halos, whereas bright objects reside in halos of > Mq. 


(vii) The halo occupation rate of LAEs and LBGs is 
very similar, except for the bright sample, where LAEs 
appear to be about four times rarer than LBGs, in broad 
agreement with the observations of Kovac et al. (20071. 


(i) We can reproduce the luminosity functions of 
LAEs and LBGs between z « 3 and z « 7 reasonably 
well. At z > 6, the model may slightly underpredict the 
abundances of sources, as it only matches data with the 
lowest observed number densities. We note that the data, 
especially for the Lya LFs, are not always homogeneous 
due to different selections and degrees of contamination, 
so it is not straightforward to constrain models very 
accurately, in particular at the highest redshifts. 

(ii) Applying selection criteria similar to those used 
in the observations, the UV LFs of LAEs predicted by the 
model are in rather good agreement with the observed 
ones, except for LAE samples selected above high EW 


(viii) More massive halos tend to host weaker 
Lya emitters on average in our model, which suggests 
that LAEs selected with high EW cuts will preferentially 
probe LAEs in low mass halos. 

In this article, we have shown that our model can 
reproduce many observed statistical quantities at z « 
3-7, using a rather simple modelling of galaxies and 
Lya radiation transfer. Some additional ingredients may 
still be necessary to accommodate the tail of the EW dis¬ 
tribution, possibly due to burstiness, and the exact scaling 
between UV luminosities and Lya equivalent width, by 
increasing the effect of dust in UV-bright objects. Nev¬ 
ertheless, within a single and coherent framework, we 
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are able to interpret reasonably well the abundances of 
LAEs and LBGs, as well as the UV LFs of LAEs and the 
Lya fraction in UV-selected samples. 

The picture emerging from this study of the Lya, 
UV, and halo properties of high redshift galaxies is con¬ 
sistent with the idea that LAEs are a subset of the LEG 
population. The apparent differences between these two 
populations would simply arise from the EW selection in 
LAE surveys, and the continuum detection limit used to 
select LBGs. In other terms, LAEs undetected in the UV 
should always be probed by deeper LBG surveys. At given 
SFR, LAEs and LBGs should have very similar UV mag¬ 
nitude, and the effect of dust on Lya photons will redis¬ 
tribute the most massive, intrinsically Lya-bright, galax¬ 
ies at fainter fluxes and lower equivalent widths, although 
these can still appear as LBGs ( |Shapley et al.|[2003[ ). As 
a consequence, LAEs that are observed are more likely to 
have lower dust extinction than LBGs. This picture seems 
consistent with the comparison study of the mid-IR prop¬ 
erties of LAEs and LBGs | |Yuma et al. ] |2010[ ), the work 
of Verhamme et al. (20081 based on stellar population 
and Lya radiative transfer modelling, the hydrodynamical 


simulations of Dayal & Ferrara (2012 I, and the results of 


Cooke 120091 and Schaerer et al. (20111 who demonstrate 
that the Lya emission properties of UV-selected galaxies 
can be statistically inferred from broad-band photometry 
and SED modelling techniques (see Dunlop 2013 for a 
review). 

Nonetheless, LAEs remain quite powerful at probing 
faint (i.e. low-mass) high redshift star-forming galaxies, 
especially strong emitters which would require very deep 
observations to be detected in UV continuum surveys. 
In addition, the observational and theoretical study of 
the Lya line profile and morphology of LAEs hold great 
potential in probing the distribution, the content, and 
the kinematics of the gas in the interstellar medium, 
but also at larger scale in the circumgalactic and in- 
tergalactic media (see e.g. [Barnes, Garel fc Kacprzakj 


20141. Increased numerical power along with larger data 


samples extending to fainter luminosities and higher 
redshift, expected from ongoing and forthcoming surveys 
(e.g. MUSE, HETDEX, KCWI or Hyper Suprime-Cam), 
will undoubtedly help refine our understanding of 
Lya emitters, and galaxy formation and evolution at 
high redshift in general. 
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